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Abstract 

We discuss a model for heavy mesons where the light quark (u or d) moves in the colour 
electric field from a heavy quark (c or b) placed in the center of the bag. We calculate energy 
spectra for pionic and photonic transitions from excited states. The transition amplitudes and 
the branching ratios between electromagnetic and pionic transitions compares favorable with 
the limited amount of known experimental data. 



1 Introduction 

Heavy quark spectroscopy is a very interesting and rewarding subject for study. The discovery of 
charmonium definitely swept away all doubts the physics community had about the existence of 
quarks as the fundamental building blocks of hadrons. 

Mesons with one heavy and one light quark is a further excellent laboratory to test our ideas 
about strong interactions. These mesons are in a way the hydrogen atoms of quark physics. As the 
mass of the heavy quark increases, its motion become gradually less and less important and the 
physical properties of the heavy- light, Qq, meson are more and more determined by the dynamics 
of the light quark. 

The discovery of the heavy quark symmetries by Isgur and Wise |^] and the creation of a 
heavy quark effective theory from QCD|| |], ^ has been extremely important for the analysis of the 
physics of heavy hadrons Q 

Ideally one would like to compute the couplings in the baryon and meson Lagrangian from QCD 
- in time this should be provided by lattice QCD calculations. In the meantime model calculations 
can be useful and one might hope that these give us some physical insight for the long distance 
behavior of the quark interaction. 



2 The model 

There are many models used in quark physics and we have chosen a variant of the M.I.T. bag model 
that was created by W. Wilcox, O. V. Maxwell and K. A. Milton (WWM), at a time when there 
were little information about exited systems made of one heavy and one light quark. 

The model is a nice theoretical laboratory, it lends itself to analytical calculations and it seems 
to give results that are not too far from experimental results. In particular it seems to work well 
for calculations of the Isgur- Wise function Q and to represents an improvement over results coming 
from the M.I.T. bag model ||. In the WMM-model the heavy quark is placed in the center of the 
bag and the light quark moves in the colour electromagnetic field set up by the heavy quark. 
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The Hamiltonian for the light quark is then 

H = H + Hi where (1) 

H = a-p + f3m + gti a V a , t a = ^ a = 1, . . . ,8 and (2) 

Hi = -gt la a ■ A a . (3) 

t a are the generators of the SU(3)c colour group. The index 1(h) refers to the light(heavy) quark. 
We have used A£ = (V a , A a ) and the usual notation a, = 7°7i and (3 = 7 , where the 7's are the 
Dirac matrices and A a the Gellman matrices. V a and A a are the colour electric potentials and vector 
fields respectively produced by the heavy quark. 

As gluon selfcouplings are neglected and the heavy quark is treated as point like the potential 
has a Coulomb like form : 

Va = ^ (4) 

Substituting this potential into equation (||) give us : 

H a = a -p + (3m+^P^. (5) 

Using the constraint that the meson is a colour singlet, that is t"t/ la = —4/3 the equation of motion 
of the light quark in the meson rest frame is then simply 

H Q = ap + 0m- ^, (6) 
r 

where £ = ia s = If-. 

The four component wave function ip(r) of the light quark (ignoring Hi) is therefore the well 
known solutions for the relativistic Coulomb problem. We shall use the notation 

where \k IS the two component spinors describing the angular part of the wavefunction. 

The energy of the confined light quark is determined by the Bogolioubov-MIT 
boundary condition (h], |ll[ which in the rest frame of the meson takes the form 
—i(f ■ ~f)ip = ip. Substituting equation (Q) into this equation and using the following property 
O • r)x£ = -X-k give us : 

f(R)+g(R) = 0, (8) 

where R is the radius of the spherical bag. The confinement of the light quark presumably originating 
from the gluonic selfcouplings is now taken care of by equation (^) and the surface conditions : 

r ■ E a = (9) 
fxB" = (10) 

The vector fields A a that are set up by the heavy quark and fulfill the boundary conditions ( px| ) 

1 m a xf m a xf 

a » = ^^^ + ^rT^ (") 

where m a is the colour magnetic moment (s) of the heavy quark : 

m a = 9 -^S. (12) 



are 
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Figure 1: Energy levels inside a bag with radius R — 5GeV 



Mh is the mass of the heavy quark and S its spin operator. Hi now takes the form : 

H I _s ! gs: 8 .( BxW i + ^, (13) 

The contribution of Hi to the energy is calculated perturbatively and the hyperfine splitting energy 
to first order in a s is 

e - = i^4^i^ F + " - ]{j + » - r *< 2 + ^ < i4) 

N is the normalization of the wavefunction, N = drr 2 (\f(r)\ 2 + \g(r)\ 2 ). Here F is the total 
angular momentum of the mesonic system and J is the light quark (total) angular momentum. 
From equation ( |l4| ) we see for Mh oo then Ej — > 0; this is the heavy quark limit. 
The mass functional for a heavy meson described in our bag model will be : 

M = M (R) = E Vo i + Ezero + m Q + E q (15) 

where E Vol = ^-BR 3 is the energy needed to create a bag in vacuum , Ezero is the zero point 
energy proportional to 1/R, tuq is the heavy quark mass and E q is the light quark energy E q = 

^jv\ + m q + E], where _E)is the hyperfine splitting energy to first order given in equation (flil). We 
will say more about the first two terms later. For now we only note that if the radii of two mesons 
with same flavour of the heavy quark is kept constant, then the mass difference between them are 
given by the following formula : 

AM = E q {nLj) - E q {n'L'j,). (16) 

This means that if the radii of two mesons do not differ too much, then the difference between the 
energy levels is directly related to the mass difference of the two mesons. 

It is of interest first to sec how the energy levels of the meson are ordered in the heavy quark 
limit when the colour electric central potential increases, these are shown in fig. [l| 

As we can see an increase in the central field from the special case where the light quark moves 
freely in the bag, reduces the mass of the P states - and even make them cross. The odd parity 
states stay however roughly half way between the ground state and the first excited 5-state. This 
will lead to some difficulties when we try to fit the spectrum of heavy meson states. 
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The breaking of the heavy quark limit is given by the spin-spin interaction and in our model by 
the term Ej, given in equation (p^). It should be noted that this term is dependent of both X/M 
and a s which determines the strength of the central four vector potential the light quark moves 
in. In this respect our model is more constraining than most models where the interquark central 
potential is unrelated to the strength of the spin-spin and spin-orbit interaction. 




Figure 2: The energy of a (massless) light quark inside a Coulomb bag with radius R = 5GeV 
and a mass of the heavy quark mt — A730MeV 



In figure g we have plotted the energy levels for the light quark in the case where M = nib for 
constant R. In this graph we have also plotted the heavy quark limit, the energy level above the 
heavy quark limit is for each pair of states where the spin of the light quark and the heavy quark 
couples to S = 1 and the level below for the case where the spins couples to S = 0. 

We see that for a given heavy quark multiplet, the induced splitting of the formerly degenerate 
states with the same angular momentum J of the light quark, but with different angular momen- 
tum F for the meson, is a highly nonlinear function of a s . Only for a s smaller than 0.2 can the 
hyperfine splitting with a reasonable approximation be taken as a linear function of a s as it is in 
the nonrelativistic quark model treatment. 

From figure || we also note that the hyperfine splitting (for finite £) increases as we go up in light 
quark excitations. This is quite opposite from the situation in the hydrogen atom and is a reflection 
of the bag models abrupt confinement. 

For the charm sector only for the two lowest states it is reasonable to calculate the hyperfine 
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splitting perturbatively, this is because of the much smaller mass of the c-quark. 



3 The mass functional for heavy mesons 

In the previous section we looked at the qualitatively features of the spectra of the heavy mesons. 
Now we will try to reproduce the quantitatively measured masses. In figure || and ^ we have showed 
the observed spectra of the D and B mesons. The mass formula for heavy mesons in our model is : 

M = ^BR* + | +m Q + y/p* + m q + ^^-^(FiF + 1) - j(j + 1) - \)1 (17) 
where T is an integral over the radial wavefunctions 
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R 



Z=xJ o dr(2 + —)f(ry g (r) (18) 

When calculating the Coulomb potential in classical electrodynamics it is common to choose 

V(r) = for r — > oo (19) 
Inspired by this we adjust the potential inside the bag to be zero at the surface of the bag [T^| : 

V(r) = -i — -^(- h (20) 
r r K 

Now the value of the potential at the bag surface will be zero, independent of the radii of the different 
mesons. Because of this transformation the light quark gets a contribution £/R to the energy. This 
we include in the light quark energy. 

The second term, C/i?, in equation ( |l7| ) is supposed to represent the zero point energy in the 
bag. When one quantizes a radiation field there will always be an infinite zero point energy term, 
but since physical quantities often are energy differences, the zero point energy falls out. However, 
when the quantization is carried out in a finite cavity, as in our model, there will be additional pieces 
of the zero point energy which depends on the size of the cavity. This is represented by the term 
C/R. The constant C can be calculated if we believe that the zero point energy give rise to the 
Casimir energy. The Casimir energy inside a perfectly uncharged spherical shell has been calculated 
by K.A.Milton, L.L.DeRaad, Jr. and J.Schwinger in ||. They obtained a value E = 0.09235/(2i?). 
To find the value in our model we simply have to multiply the value by eight, because there are 
eight gluonic radiation fields : 

C-0.37 (21) 
We have determined the masses by minimizing the mass functional by the relation : 

It turns out that the integral in equation (|l^) goes as 

tiiqR 2 ^ ^ 

Since p q ~ 1/R it is clear that when R — ► 0, T will dominate. If the factor in front of the hypcrfinc 
splitting is negative, equation ( p"7| ) will for some choice of the parameters have no finite minima. 
When the radius become small we can not neglect the repulsion of the heavy quark and the model 
becomes meaningless. This is a well known problem and the usual way of dealing with this is to 
argue that the divergence will disappear when we calculate higher orders correction. The procedure 
then is to minimize the energy with respect to R before adding the hyperfme term Jl4[ . In our model 
states with same J and L will then have equal radii. 
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4 Heavy meson masses 



First of all we have to determine the parameters in the model. There are four parameters the 
strong coupling constant (a s ), the bag constant(-B), the heavy meson mass(mg) and the light quark 
mass(m 9 ). There are of course many ways of determining the parameters, we have chosen too look 
at the B mesons. We have chosen to determine a s ,B,m,b from the observed masses of the B. B* 
and B' mesons, and since the u and d quark have a small mass we have assigned to them a (rather 
unimportant) mass of lOMeV. The calculated results are shown in table |l| 



Experimental 
Results 


Theoretical 
Results 




Mass [Me V 




M&ss[MeV] 




State 


RadiusIGeV^ 1 ] 




B 


5279.4 ±2.2 [ 


15 




5279 


o- 


lS'l/2 


3.94 





B* 


5324.8 ± 1.8 [ 


15 




5325 


1- 


lS'l/2 


3.94 


1.52e, + 0.203eg 


B 


? 


5592 


0+ 


2A/2 


4.50 





B! 


? 


5671 


1+ 


2Pi/a 


4.50 


0.625e g + 0.203e Q 


B! 


5725 


ie 


1 


5623 


1+ 


2^3/2 


4.47 


1.93e 9 + 0.101e Q 


B 2 


5737 


1C 


1 


5637 


2+ 


2^3/2 


4.47 


2.33e 9 + 0.203e Q 


B' 


5859 


17 




5859 


CP 


2S'i/2 


4.90 





B'* 


? 


5967 


1- 


2<S 1 /2 


4.90 


0.597e 9 + 0.203e Q 



Table 1: The parameters are B 1 / 4 = 161MeV,£ = 0.538 (a s = 0.404), m u = m d = lOMeV and 
nib = A627MeV.e q = eQ = —1/3 for d and b quarks and e q = cq = 2/3 for u and c quarks. 



In table [l] there are only given uncertainties for the IS1/2 states. This is because the value of the 
masses for the 2P 3 / 2 and the lowest 2S\/2 states was found in |l6| and |l7j as a fit to the experimental 
data and it is not quite clear how to put errors to these numbers. 

Unfortunately we see that the calculated masses for the B\ and B2 mesons are almost exactly 
lOOMeV below the experimental values. However, we see that the splitting between the two P 3 / 2 
mesons is quite well described. One may wonder if there are another set of parameters which will 
give us a good fit to all the experimental masses listed in table |l|. We believe that this is not so. The 
parameter in the system which has most influence on the spectra is £. In figure |l| we plotted the 
energy levels for £ between and 0.9. The main problem is that the P states does not move much 
in direction of the 2S , 1 / 2 state when £ increase. The boundary condition (g) clearly gives a wrong 
level splitting. A hope for the model would be that the state we used as B' is not the right one. If 
we used the B\ and B2 states (together with B and B*) as a normalization the resulting B' would 
have a mass around 6020MeV r . 

Now we have the value of the strong coupling, the bag constant and the b quark mass. In order 
to calculate the D and B s meson masses we have to find the c and s quark masses. This has been 
done by demanding that : 

m( D Ex P .) - m(D T eor) ~ m(D* Theor ) - m(D* Exp ) (24) 
and the s-quark mass by demanding that: 

m(B sE xp.) - m(B sT heor.) ~ m ( B * sTheor ) - m{B* Exp ) (25) 
The results are listed in table ||. 



5 Decays 

We have looked at electromagnetic and pionic transitions between mesons listed in table [j] and] 
The pionic transitions are calculated using the surface coupling version of the chiral bag model [pi 
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Experimental 
Results 


Theoretical 
Results 




Mass[MeV] 


Mass[AfeV] 


F 


State 


Radius^e^" 1 ] 




D 


1866.9 ±0.7 


15 




1854 


0" 


lS'l/2 


3.94 





D* 


2008.4 ±0.7 


15 




2022 


1- 


lS'l/2 


3.94 


1.52e g + 0.725e Q 


B s 


5369.6 ±2.4 


15 




5363 


o- 


l£l/2 


3.91 





b; 


5416.3 ±3.3 


15 




5424 


1" 


lS'l/2 


3.91 


1.52e 9 ±0.203eQ 




? 


5667 


0+ 


2Pl/2 


4.47 





Bi s 


? 


5737 


1+ 


2Pl/2 


4.47 


0.619e 9 + 0.203e Q 


Bi. 


5874 | 


16 




5718 


1+ 


2^3/2 


4.45 


1.93e 9 + 0.101e Q 


B 2s 


5886 | 


16 




5732 


2+ 


2^3/2 


4.45 


2.31e 9 ±0.203e Q 


B' s 


? 


5887 


0" 


2S'i/2 


4.87 





B's 


? 


6008 


1- 


2S'i/2 


4.87 


0.593e 9 + 0.203e Q 



Table 2: The parameters are B 1 / 4 = 16LMeF,£ = 0.538 (a s = 0.404), m b = A627MeV, m c = 
129AMeV and m s = 231MeV. e q = eQ = —1/3 for d and b quarks and e q = eQ = 2/3 for u and c 
quarks. 



In this model a pion field carries the axial current outside the bag produced by the light (u or d) 
quark inside the bag. This makes the model chirally symmetric for massless u and d quarks and the 
interaction between the bag and the pion field is given by : 



Ant — 



2/. 



lpJ 5 T • 4>ipA s , 



(26) 



Tj are the Pauli isospin matrices and cf> an isovector representing the pion field. A s is a covariant 
surface delta-function. 

The calculation of transitions involving pions is then straight forward, some expressions for pionic 
transitions are listed below : 



1"0" 



+ 



0~0~ 



1"0" 



1"0" 



1"0" 



o _ o- 



^p(R) (jr^k) + Y^k) + r°(fc)) a 



v fi 

r{B'^B**) = ±%\P(R)C v \ 



T(B> — B n) = ^\S(R)C n \ 2 



Vfi 

r(J 

Vfi 



r(B 2 ^Bir) = ^jt\D(R)C„\ 



T(B 2 — B*n) = ^-^\D(R)C n \ 2 



Vfi 

r(5x(P 3/2 ) — » B*ir) 



tjrD(R) (V6Yi(k) - V3Y 2 \k) + y 2 °(fe) a 



J_ k 



D(R)C„\ 



Vfi = j;^v^S(R)Y (k)C„ 
r(Bi(P 1/2 ) — » B*ir) = -^Jl\S(R)C* 



T(B — Bw) = ^\S(R)C^ 2 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 
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jp a _^ JPH^ 


a — ► (3tt 


Teor.[Me^] 


9a(3-K 


0" > 1-7T 


f>> k R*-7T + 


84 1 


49 8 

4Z .O 


R' . R*-7T U 


49 9 


OU.l 


0- > 0+7T 


R' > Rn7T + 


4^ 9 

4:0 .O 


91 n 


ft' > R n7r 


99 7 


1 fi 9 


2+ > l-TT 


R v R*7T + 


n 17 
y . i / 


UU.O 




4 fi7 


1Q 4 
oy .4 


2+ > Q-n 


R+ > ru„+ 

2 


Q 74 


41 ^ 


Rl) y RU„U 

-D2 ^ i J /l 


4 Q4 


10 7 
OU. 1 


R+ > R+^U 

2 


4 Q8 


10 7 

OU. 1 


Rl) > R+7T- 


Q 89 


41 ^ 

40.0 


l + (P 3 /2) — » l-TT 


, R*-Tr+ 


1 1 1 


79 8 


, R*7T U 


fi 71 
u. t 1 


Ol . 1 


l+fPi/o) ► 1 - 7T 


Bl > S*7T+ 


92.7 


27.1 


Si > £*7T U 


46.1 


19.0 


0+ — ► O^vr 


B + — > B U 7T+ 


93.7 


28.5 


5< J — > B U 7T U 


46.8 


20.1 


5 + — ► S+vr" 


46.9 


20.1 


B{j — > B+7T- 


93.7 


28.5 


1- — ► O^TT 


£>*+ ► L> u 7r+ 


6.01 10" 2 


17.1 


D*+ — > D+7T U 


2.72 10" 2 


12.1 


£* u — ► £ U 7r u 


3.88 10" 2 


12.1 



Table 3: Pion decays 



tro - 



T(D* — 



^-P{R)Yl{k)C v 

D *) = sfl p ( fl ) c >l 2 



(34) 



In these formulae 



S(R) =R 2 (r (R)g a (R) + g* (R)f a (R)) j (kR) (S-wave.) 
P(R) =R 2 (f* (R)g a (R) + g* (R)f a (R)) n (kR) (P-wave.) 
D(R) =R 2 (f;(R)g a (R) + g* (R)f a (R)) j 2 {kR) (D-wave.) 



The index a(f3) refers to the initial(final) meson. 



(35) 
(36) 
(37) 



for 7T 

for 7T°. 



(38) 



The above transition rates ((27j)-(|34D have been numerically calculated and are listed in table ^. 

In addition to the partial widths listed in table || we have also calculated the coupling constants 
(to the right in table 0) , by using the following definition : 



9a(3-K 



r(a — ► /37r)247rAf2 



k 2L+l 



(39) 



k is the pion momenta and M a the mass of the decaying particle. The dimension of the couplings 
goes as [GeV]~ L+1 , where L is the relative angular momentum between the decay products. Only 
L = 1 transitions such as 1~ — > 0~7r are then dimension less by the definition (p9|). 
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a — > fiir 


Tcor.[fceU] 


0" — > 2+ + 1- 


• D * -°27 


2.62 10~' 5 


R** + > R+^ 


1 05 1 fi~ 2 


0" — > 1+ + 1- 




1 77 1 O 1 


R** + > R+^ 


7.07 10 1 


o- — > 1" + 1" 


£**u > _B* U 7 


1.20 10 _1 


5**+ > _R*+7 


7.89 10 _i 


o- ► 0~ + 1~ 


► _R U 7 


o 


5**+ > 


o 


2+ — > 1+ + 1~ 


5<> ► 


6.26 10~' 2 


Bn ► Bn + "7 

"2 1 / 


2.50 10~ a 


2+ > 1~ + 1~ 


RJJ > R*<% 


6.76 10 1 


— ► 


2.70 10 2 


2+ ► 0" + 1~ 


#y > s'S 

±J 2 ^ 1 


4.08 




1.63 10 1 


2+ — ► 1+ + 1~ 


B" — >B ,7 


3.38 10~ 4 


2+ — > 1- + 1" 


"s2 "s / 


5.86 10 1 


2+ — ► 0" + 1- 


"^2 u « / 


4.62 


1+ — > 1- + 1" 


R<> > R*U^ 


9 7^ i n 1 

Li, 1 O 1U 


R+ > R* + ^ 


1 HQ 1 C\' z 


1+ — > 0~ + 1" 


B<> — > S° 7 

1 / 


4.10 10 1 


B+ y S+7 


1.64 10 2 


1+ — ► 1~ + 1" 




2.56 10 1 


1+ — > 0" + 1- 


^ — 3?7 


3.46 10 1 


l- — > o- + 1- 


S* u — ► B u 7 


6.41 10~ 2 


B*+ y B+ 1 


2.72 lO" 1 


Bf — ► B^'7 


5.10 10~ 2 


D *u — , D » 1 


7.18 


D*+ — y D+ 1 


1.73 



Table 4: Photon decays 



It is also possible to calculate the coupling g B " g T , the B* emits a virtual pion. This coupling has 
been calculated in the rest system of the heavy meson at zero recoil, the result is : 

9b*b^ = \[\t\ r2 imR)9a(R) + g* (R)f a (R)) \M B * (40) 

9b*b^> = \j\jJR 2 {ta(R)g*(R) + 9}{R)UR)) \m b .. (41) 

Using the wavefunctions for the B* and B meson give us 

9B*B,+ = 45.6 (42) 
g B , B7[ o = 32.2. (43) 

The coupling of photons to the mesonic states are done straight forward by using the interaction 
Lagrangian : 

Ant = e ? e^7 • Aip (44) 
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Some expressions for electromagnetic transitions are listed below 



2+1- r= fae 2 q u; k \fdxj 2 (F + G)\ 



4 1+1" 
* 1 _ 1" 

2+ — > i+r 
2+ — > i~r 

2+ — ► _ 1" 

i+ — > l-r 

i+ — > 0"1" 
1" — > o~v 

where we have defined : 



>e 2 q LJ k ( \ \Jdxj 2 (F + G)\ 2 

\Jdx(Fj 2 - G J0 )\ 2 +§ \Jdx(Gj 2 + Gj ) 



faelcJklfdxjxiF + G) 



2 
375 



ae 2 q u> k (f IfdxhiF + G)] 



-52 |/da:(Gj 3 + Gh)\ + 2 |/da;(Gj3 - Fj x )\ 
-27 j/ da;(Fj 3 - Gjt)\ 2 + 27 |/ dx(Fj 3 + F n ) 
fdxj 3 (F + G) 



3501 
7 



r = §ae*u; fe \fdxj 2 (F + G)\ 2 + |/£fe(Gia - Fj )\ 2 
+l\Jdx(Fj 2 + F JO )\ 2 ) 

?= |ae^ fc |/dxj 2 (^ + G)| 2 

r = |a e ^ fc (l3 |J dxj 2 (F + G)\ 2 + § | f dx(Fj 2 + Fj )\ 2 
+2\Jdx(Gj 2 -Fj )\ 2 ) 

r = fae|w fc (-i |/rf^ 2 (F + G)| 2 + ± |/ <fe(Fj 2 + Fj )\ 



+ \J dx(Gj 2 -Fj ) 
\a.e\w k |/ dxji(F + G)| 2 

^ = x 2 fp{x)*g a {x) 
G = x 2 gp(x)* f a (x) 
ji = ji(oj k x) 



(45) 

(46) 

(47) 
(48) 

(49) 

(50) 
(51) 
(52) 

(53) 
(54) 



(55) 
(56) 
(57) 



The numerical values for the above expressions are shown in table [|. The transitions in table 
^ are very suppressed in comparison with those listed in table |j| This of course is expected, the 
smaller phase space for pion decays is compensated by the much stronger pion coupling relative to 
the electromagnetic coupling. 



5.1 Comparison with theoretical and experimental results 

We have calculated a lot of partial widths for different particles, listed in table || and f|. To day very 
little is known on the experimental front, but there are a lot of theoretical predictions. So we will 
compare our results with the known experimental and some of the theoretical results. As we shall 
see there are no conflict between our predictions and the experimental information. 

In table § we have listed some theoretical and experimental results, on the experimental limits 
we have assumed that the the width of D*° have the same upper limit as the width of D*^. It may 
not be clear from table |[ but the theoretical predictions vary a lot. In |22| there is a summary of 
theoretical estimates. For the particular decay D* + — ► D°ir + that determines the coupling constant 
9d*Dtt, the predicted rates vary from 10fceV^(QCD sum rules) to more than 100fceV(quark model + 
chiral HQET). We obtained T(D*+ — ► D q tt+) = 60.lkeV, close to the value (61 - 78keV) coming 
from P.Cho and H.Georgi who makes calculations with chiral HQET. The value of the coupling 
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Theoretical 


Experimental 


HGIfMeV") 


l2ClffeeV") 


\2l\(keV) 


our work 


\l&{keV) 


T{B 2 — ► B*ir) 


11 






l3.8MeV 




T{B 2 — ► Btt) 


10 






U.7MeV 


— 


T{B 1 — > B*tt) 


14 






20.0MeV 




r(p*+ - -> p 7 ) 






0.38 ±0.06 


0.272keV 




r(s* u — ► p u 7 ) 






0.13 ±0.03 


6.41 lQ-' 2 keV 




r(p* — > p s7 ) 






0.22 ±0.04 


5.10 10-' 2 keV 




r(Z3*+ — ► D u tt+) 




69.1 




m.lkeV 


< 91 


r(D*+ — ► d+tt u ) 




32.1 




27.2keV 


< 44 


r(£* l) — ^ d u tt u ) 




46.0 




38.8keV 


< 85 


r(D*+ — > l>+ 7 ) 




0.919 


0.23 ±0.1 


l.72keV 


< 4.1 


r(z?* u — -> D u 7 ) 




23.5 


12.9 ±2 


T.lSkeV 


< 54 



Table 5: Comparison with theoretical and experimental results 



9b*Bit+ vary from gB*Bn+ = 15 ± 4(QCD sum rules) to gB*Bn+ — 64(quark model±chiral HQET), 
we found gs*Bir+ = 45.6. 

We have calculated most but not all decay modes for the excited states. The irir modes are 
missing. As we believe that these modes are less important than the emission of single pions in the 
decays we still can give approximate values of the decay widths, these are : 



251/2 
2^3/2 
2Pl/2 



T(B') ~ 195MeV (58) 

T(B 2 ) ~ 29MeV r(B 1 ) ~ 20MeV (59) 

r(S x ) ~ !39MeV T{B ) ~ UlMeV (60) 

T(D* + ) ~ 89keV T{D*°) ~ 46keV (61) 



The P1/2 states are naturally much wider than the P3/2 states because they decay only trough an 
S-wave, whereas the P 3 / 2 states decay through a D-wave. The full width of the P states indicated 
by a preliminary experiment Jji) are T(B(P 3/2 )) ^ 20MeV and L(P(P 1/2 )) ~ 150MeV. This is, 
as we see, in good agreement with our results. Since the Py 2 states are so broad, they are very 
hard to reconstruct from the experimental results, and so far there have not been any really precise 
measurement of their masses. 

The CLEO report |^5| contains best measurement of the D* + branching fractions, a large im- 
provement of what can be found in the Particle Data Book |l5j . Since we have calculated the width 
of the D* + meson, equation (^lj), it is easy to calculate the branching ratios. The results are shown 
in table || together with the CLEO results. Our results are clearly in good agreement with the ex- 





CLEO 


our calculations 


B r (D*+ — > L>+ 7 ) 


(1.68 ±0.51)% 


1.94% 


B r {D*+ — > D+ir {> ) 


(30.73 ±0.63)% 


30.55% 


B r {D*+ — ► D+TT+) 


(67.59 ±0.70)% 


67.51% 



Table 6: Comparison with CLEO data and our predictions 



perimental data We recognize however that branching ratios are one thing, particular decay widths 
another. We get the correct ratio between pionic and electromagnetic decays. As the coupling of the 
electromagnetic field to the quarks is simple, we naturally have some confidence in the calculated 
electromagnetic transitions rates. Therefore we believe that our calculated pionic rates cannot be 
too far off from what will be measured in the future. 
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Figure 3: D meson spectra and transition lines. 
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Figure 4: B meson spectra and transition lines. 
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